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Abstract 

The production of the Standard Model Higgs boson in the four- 
fermion reaction e^e~ — > fi'^fi^bb is studied. The complete tree-level 
matrix element, including signal and backgrounds in the standard elec- 
troweak theory, is computed and initial state radiation is taken into 
account in the leading-log approximation. A Monte Carlo event gen- 
erator has been built and numerical results for some distributions of 
experimental interest for the search of the Higgs particle at future 
electron-positron colliders are shown, compared with those existing in 
the literature and commented. 
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After the experimental confirmations of the standard theory of the elec- 
troweak interaction provided by the precision tests at LEP/SLC and by 
the evidence of top quark production at the TEVATRON the search 
for the Standard Model Higgs boson becomes a primary task of the exper- 
iments planned at future electron-positron and hadron colliders. For the 
case of a light or intermediate Higgs particle {Mh < 2Mw), the Bjorken 
bremsstrahlung process e~^e~ — * HZ is known to be the main production 
mechanism at e~^e~ colliders up to 0.5 TeV |^. Detailed studies about the 
production of the Higgs boson in e+e~ collisions are given in Refs. @, |]. 
From these investigations it emerges that the specific four-fermion reaction 

e+e" ^ HZ ^ fi+fi-bb (1) 

provides the cleanest event signature coupled with "small" backgrounds. 

A complete tree-level calculation of the four-fermion process (|ID, includ- 
ing all possible diagrams and their interferences, has been performed in two 
recent papers 0. On the one hand, in Ref. a Monte Carlo approach 
has been adopted but the effect of the initial state QED corrections, large 
and hence necessary for the data analysis, has been neglected; on the other 
hand, initial state radiation has been included in 0, but according to a 
semi-analytical procedure which is limited from the point of view of event 
simulation purposes. The aim of this letter is to present a Monte Carlo event 
generator for the process e~^e~ ^^^~bb based on the calculation of the 
full set of diagrams contributing in the standard electroweak theory to the 
assigned four-fermion final state and on the inclusion of initial state QED 
corrections in the leading-log approximation. We will only sketch the basic 
theoretical ingredients, putting more emphasis on the phenomenological as- 
pects of the study and on the potentials of a "specialized" event generator 
for future searches of the Higgs boson at LEP2 and NLC energies. 

The kinematics of the 2—^4 process at the Born level has been treated 
according to the general recipe for a four-fermion final state described in . 
Essentially, seven independent dimensionless quantities are introduced and 
chosen, together with a trivial overall azimuthal angle, as independent vari- 
ables from which the full four-momenta of the outgoing particles are recon- 
structed in the laboratory frame by explicit solution of the kinematics. 

The inclusion of initial state radiation is simply achieved by dressing the 
incoming leptons with QED structure functions 0. They allow to include, 
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by virtue of factorization theorems, the long-distance universal contributions 
in soft and/or collinear approximation. 

Taking into account the initial state radiation, the total cross section can 
be expressed as a nine-fold convolution of the form 



where d[PS] denotes the volume element in the 7-dimensional phase space 
connected to the independent variables of the hard scattering reaction in the 
cm. frame. The scale of the hard scattering cross section is the reduced 
cm. energy s = X1X2S which is the main effect of the initial state radiation 
together with the boost of the cm. frame with respect to the laboratory 
one, which is also taken into account by properly rescaling the independent 
variables of the hard kernel cross section in the cm. frame. 

In the kernel cross section the complete set of tree-level diagrams con- 
tributing in the Standard Model to e~^e~ — >• fi~^fi~bb is included. Twenty-five 
diagrams have to be computed and they can be classified as follows: 

• the signal diagram e~^e~ — > HZ followed by the decays H ^ bb and 

• the eight t— channel diagrams relative to ZZ, 77, Z7 exchange; 

• the sixteen s— channel diagrams corresponding to the bremsstrahlung 
of a neutral vector boson from the fermionic final states (radiative pro- 
cesses) . 



Each amplitude has been computed by means of SCHOONSCHIP |10| within 



the helicity amplitude formalism |TT[ following the general strategy described 
in 1^. This allows to express the kernel cross sections essentially as the sum 
of products of resonating propagators with functions of the dimensionless 
variables originating from the calculation of the traces of the 7 matrices' 
strings. Mass corrections have been neglected in the kinematics and in the 
calculation of the amplitudes. These contributions introduce terms of the 
order of ml/s which at LEP2 energies amount to about 0.1% and hence they 
can be neglected as long as the ratio signal / background is larger than some 
per cent. 
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To take care of the peaking behaviour of the integrand the Monte Carlo 
event generation is treated according to the importance samphng techni- 
que |T^. The latter is applied, by means of proper changes of variables, to 
the exponentiated infrared sensitive part of the structure functions, to the 
different propagators' "singularities" (particularly to the Higgs propagator 
which resembles for its very narrow width a delta-like function) and finally 
to the jacobian peaks appearing in the phase space factor. In particular, the 
weights Wi and W = Y^iWi are introduced to obtain a proper samphng of the 
double propagator "singularities" of the signal (ZH) and the backgrounds 
{ZZ, Z'y and 77 in all possible combinations). 

Let us show and comment some results obtained with our event generator 
based on the random number generator RANLUX |T^. The input parameters 
used are: Mz = 91.1887 GeV, Tz = 2.4974 GeV, Mw = 80.22 GeV, = 
4.7 GeV WM. For the Higgs boson width the tree- level expression is adopted. 



In Fig. la the total cross section in Born approximation is compared with 
the initial state QED-corrected one as a function of the cm. energy and for 
three different values of the Higgs mass. The dashed, dash-dotted and closely 
dotted lines represent the lowest-order results for Mh = 65, 100, 140 GeV, 
respectively; the three solid lines refer to the above Higgs mass values and 
include the effect of the leading-log initial state QED corrections. As a ref- 
erence curve, the QED corrected total cross section obtained by excluding 
the Higgs boson contribution is also shown (sparsely dotted line). Cuts on 
the invariant masses of the and bb pairs of M^+^- > 12 GeV and 

Mf^b ^ 12 GeV (typically used in other numerical studies |^) are apphed 
in order to reduce a large fraction of the 77 background. The values ob- 
tained for the cross sections with our event generator have been compared 
in detail (i.e. for different Higgs masses and cm. energies) with the inde- 
pendent results of Refs. 0, |^ and found to be in good agreement. It should 
be noted that in the LEP2 energy range the size of the total cross section 
of the fi~^fi~bb reaction can give by itself evidence of Higgs boson production 
in the "light mass" interval 65 GeV < Mh < 100 GeV (see also Fig. lb), 
whereas for higher Higgs mass values a cm. energy larger than 200 GeV 
is necessarily required to distinguish the signal from the background. This 
feature, already pointed out in Ref. |0] for the lowest-order cross section, is 
not largely affected by initial state bremsstrahlung which essentially lowers 
the peak cross section of about 10% and originates a moderate radiative tail 
at high energies {^/s > 300 GeV ), as generally occurs for other "double 
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resonant" processes {WW, ZZ) above threshold. Moreover, the two-peak 
shape for Mh = 65 GeV and Mh = 140 GeV is made more flat by effect 
of the initial state radiation, which however does not cancel the two-bump 
structure. Their appearing is given, with very good approximation, by the 
incoherent sum of the signal and backgrounds diagrams. In agreement with 
Ref . , we checked that the interference terms between diagrams with differ- 
ent intermediate bosons and also s— channel radiative processes give a small 
contribution. 

Figure lb shows the ratio cTs/o"/ as a function of the Higgs mass and 
for three different energies of interest at LEP2. cr^ denotes the signal cross 
section and o"/ stands for the full four-fermion cross section including signal 
and backgrounds. The effect of the initial state radiation is included in this 
plot. As already singled out by Fig. la, for the case of a light Higgs mass 
ranging between, say 65 GeV and 80 GeV, the ratio cr^/o"/ turns out to be 
more favourable for "low" centre of mass energies; on the contrary "high" 
centre of mass energies maintain the ratio more flat and sensitive to a Higgs 
mass up to, say, 100-110 GeV. 

In Figs. 2-5 we show some realistic distributions of experimental interest 
for the Higgs search at LEP2. These results have been obtained processing 
the n-tuple created by the event generator at = 205 GeV. The full set 
of tree-level diagrams and the contribution of the initial state radiation have 
been taken into account in the simulation. A sample of 10^ four-fermion 
events has been generated in order to study with small statistical error the 
behaviour of the distributions when varying the Higgs mass. Although differ- 
ent cross sections correspond to different Higgs mass values, the same number 
of events has been considered in order to study the pattern of evolution of 
the shapes while varying the Higgs mass. The hh invariant mass distribution 
(Fig. 2), the 6— quark scattering angle distribution (Fig. 3) and the hh rel- 
ative angle distribution (Fig. 4) turn out to be useful observables to detect 
a Higgs signal and, if found, extract information on the quantum numbers 
assignment. From Fig. 2 it can be seen that the hh system invariant mass 
distribution exhibits two clear spikes positioned at Mh and M^, as long as 
Mh < say 110 GeV, the spike at Mh resembling a delta-like function. The 
rise at low invariant masses is due to 7 ^ 66 in 77 and 7Z diagrams. From 
Fig. 3 it can be seen that the 6— quark angular distribution shows, for the case 
of a light mass Higgs where the signal dominates the fi'^fi~hh cross section, a 
clearly isotropic spin-zero behaviour gradually disappearing as increasing the 
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Higgs mass. The bb relative angle distribution represented in Fig. 4 shows 
a two-peak structure. The peak positioned at lower cosi^^b reflects the pres- 
ence of ZZ diagrams, while the (much smaller) one at cos ■j?^^ ~ 1 comes from 
7 — > 66 in the 77 and 7Z contributions. The peak at intermediate cosi?;,^, 
well visible in the case of a light mass Higgs, corresponds to the contribution 
of the signal diagram. The different position of the Higgs peak in the two- 
bump shape is a consequence of the different Lorentz boost acting on the 
Higgs boson depending on its mass value. The residual bump at intermedi- 
ate cos'diii for Mh — 140 GeV comes from Z ^ bb in the 7Z contributions. 
Figure 5 shows two typical photon energy distributions (for Mh = 65 GeV 
and 140 GeV, respectively) where the peak at very low energies is due to 
dominating soft multiphoton emission. 

To summarize, we presented the results of a Monte Carlo event generator 
for the four-fermion reaction e~^e~ n~^iJ,~bb, including all the Standard 
Model diagrams and the effect of initial state radiation in the leading-log 
approximation. 
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Figure Captions 



Figure 1. (a) The total cross section of e+e~ — > fjL'^ii'bb without (dashed, 
dash-dotted and closely dotted line) and with (solid lines) initial state leading- 
log QED corrections as a function of the cm. energy and for different Higgs 
masses. The sparsely dotted line is the QED corrected total cross section 
obtained excluding the Higgs signal, (b) The ratio o's/'^f as a function of the 
Higgs mass in the energy range of LEP2. (jg is the Higgs signal cross section, 
(7/ the full cross section including signal and backgrounds. QED corrections 
are included. 

Figure 2. The bb invariant mass distribution for four different Higgs masses 
at ^/s — 205 GeV (log. scale). Initial state radiation is included. 

Figure 3. The 6— quark scattering angle distribution for four different Higgs 
masses at — 205 GeV. Initial state radiation is included. 

Figure 4. The bb relative angle distribution for four different Higgs masses 
at y/s — 205 GeV. Initial state radiation is included. 

Figure 5. The photon energy distribution for two Higgs mass values at 
205 GeV (log. scale). 
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